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Abstract 

The type la supernovae observational data is one of the most important in observational cosmology 
nowadays. Here we present the first public version of BETOCS (BayEsian Tools for Observational 
Cosmology using SNe la), which is a powerful and high productivity tool aimed to help the theoretical 
physicist community investigate cosmological models using type la supernovae (SNe la) observational 
data. BETOCS is applied to the generalized Chaplygin gas model (GCGM), traditional Chaplygin gas 
model (CGM) and ACDM, ranging from 5 to 3 free parameters, respectively. The "gold sample" of 
157 supernovae data is used. It is shown that the Chaplygin gas scenario is viable (in most cases the 
ACDM is disfavoured) and the quartessence scenario (that unifies the description for dark matter and 
dark energy) is favoured. The Hubble parameter (Ho) is important and should not be fixed and it can 
be estimated or marginalized with or without the Hubble Space Telescope prior. 
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1 Introduction 

The type la supernovae (SNe la) observational data has forced us to discard or change the majority of 
theoretical cosmological models supposed to be correct until the second half of the last decade Q][2]. The 
crossing of the SNe la statistics with other observational data, like the anisotropy of the cosmic microwave 
background radiation (CMBR) 3], gravitational lensing the X-ray gas mass fraction of galaxy clusters 

[S] , etc, leads to a scenario where the matter content of the Universe is described by an unclustered component 
of negative pressure, the dark energy, and a clustered component of zero pressure, the cold dark matter. 

There are many candidate for dark energy, the most natural seems to be the cosmological constant 0, 
since it can be connected with the vacuum energy in quantum field theory [5] , but the small value resulting 
from observations for the energy density of the cosmological constant term yields a discrepancy of about 120 
orders of magnitude with the theoretically predicted value ^J. Among many other possibilities, for example 
there is the quintessence model with scalar fields |1(JI 111) . 

Here we will focus on the Chaplygin gas models (CGM) |121 1131 IT4l I15| . It is based on a string inspired 
configuration that leads to a specific equation of state where pressure is negative and varies with the inverse 
of the density [IB]. This model has been generalized, giving birth to the generalized Chaplygin gas model 
(GCGM), where now the pressure varies with a power of the inverse of the density |15|. These proposals have 
many advantages, among which we can quote the following: in spite of presenting a negative pressure, the 
sound velocity is positive, what assures stability |17j : these models can unify the description of dark energy 
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and dark matter, since the fluid can clusters at small scale, remaining a smooth component at large scales 
|15| : the CGM has an interesting connection with string theory JHj- Some criticisms have been addressed to 
the GCGM (CGM) mainly connected with its features related to the power spectrum for the agglomerated 
matter |18j . However, in our opinion, this specific criticism is not conclusive, since the introduction of baryons 
may alleviate the objections presented against the cosmological scenarios based on the GCGM (GCM) |T5). 

The GCGM (Generalized Chaplygin Gas Model) is defined as a perfect fluid with an equation of state 
given by 



where A and a are constants. When a = 1 we re-obtain the equation of state for the CGM (Chaplygin 
Gas Model), the traditional Chaplygin gas model. See refs. |2()l 1211 l2"2~] for more detailed definitions of the 
Chaplygin gas models used in the present work. 

All free parameters for each model are considered. In the case of the GCGM there are five free parameter: 
the Hubble constant Hq; the equation of state parameter a; the "sound velocity" related parameter A; the 
density parameter for the pressurclcss matter Q, m o', the density parameter for the Chaplygin gas Q c q (or 
alternatively the density parameter for the curvature density of the Universe fifco). For the CGM, the 
number of parameters reduce to four, since a = 1. We also consider the ACDM, where the number of 
parameters reduce to three: H , f2 m o and Q c0 (alternatively, £lko)- 

One important point is how to perform this statistical analysis: the final conclusions may, in some cases, 
depend on the statistical framework (Bayesian, frequentist, etc.), as well as on the parameters that are 
allowed to be free, and how these parameters are constrained (through a joint probability for two parameters, 
minimizing the error function or through a marginalization of all parameters excepted one, etc.). In some 
cases, the different procedures adopted may lead to quite different conclusions on the best value for a given 
set of parameters. The choice of the observational data sample may of course be important as well. 

The present work is intended to : 

1. Announce the first public version of BETOCS (BayEsian Tools for Observational Cosmology using 
SNe la [231 )i which is a powerful and high productivity tool aimed to help the theoretical physicist 
community investigate cosmological models using type la supernovae (SNe la) observational data. BE- 
TOCS is a freeware and open source tool written in the Mathematica |24| language. The Mathematica 
notebooks of BETOCS contain documentation, source code and practical examples with textual and 
graphical outputs. 

2. Emphasize that fixing H Q is not acceptable, yielding arbitrary parameter estimations and usually bad 
best-fittings. On the other hand, the HST (Hubble Space Telescope) prior for H implies minor 
effects on all best-fittings and parameter estimations (as shown here comparing the tables with flat and 
HST priors for Hq), so its use is just a matter of choice. Nevertheless, it is recommended to use the 
HST prior as it is standard when using other observational cosmological data (X-ray gas mass fraction, 
etc). 

3. Continue the work of refs. |2(JI 1211 122] and show a very complete view and analysis of the Chaplygin 
gas models (generalized and traditional) using SNe la. The present work specially fixes an error on ref. 
|2"U] when calculating the best-fittings and parameter estimations for non-flat Universes ( flko ^ 0), 
due to an earlier version of BETOCS with wrong optimization code. B ut the side-effects of this error 
were not critical, just worsening some positive features. 

The best-fitting parameters and the parameter estimations were calculated for each case of GCGM, 
CGM or ACDM using a corresponding BETOCS (BayEsian Tools for Observational Cosmology using 
SNe la [22]) notebook (written in Mathematica [2J] language) containing the : definitions of the theoretical 
cosmological model, reading of the observational SNe la data, \ 2 definition, Bayesian tools library, calculation 
of the Bayesian PDF, global maximization of PDF, PDF visualization in 3 dimensions (if available), PDF 
visualization and analysis in 2 dimensions and finally PDF visualization in 1 dimension with parameter 
estimation. This article does not include all the graphics and analyses of the BETOCS notebooks, but some 
of them are available on the Internet site of the BETOCS project |23| . 
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This paper is organized as follows. In section we detail the best-fitting analysis using BETOCS, with 
many results presented in tables. SectionOlexplams how the parameter estimations are made using BETOCS, 
such that a detailed Bayesian analysis is performed for each independent and dependent parameter, with 
the results shown in many tables and figures. The conclusions are discussed in section 0J 



2 Best-fitting parameters using BETOCS 

In order to compare the theoretical results with the observational data, the first step in this sense is to 
compute the quality of the fitting through the least squared fitting quantity \ 2 - In the case of flat priors for 
all independent parameters of the theoretical cosmological model, we get 

X 2 = E (/i "%^ )2 • (2) 

In this expression, fig i is the distance moduli observationaly measured for each supernova of the 157 gold 
SNe la dataset Mo % ^ s the value calculated through the theoretical cosmological model, er^ o i is the 
measurement error and includes the dispersion in the distance modulus due to the dispersion in galaxy 
redshift due to peculiar velocities, following ref. (20] • It is useful to define xt '■ X 2 divided by the number of 
degrees of freedom of the observational data, i.e., the number of SNe la, here 157. 

As we also want to compare the fitting and estimation of the parameters without priors and with the 
HST (Hubble Space Telescope) prior 25 for Hq, then the \ 2 used for the calculations with the HST prior 
is simply 

2 _ it^A ~ _l_ (-^0 - 72 ) 2 fo\ 
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Table 1: The best- fitting parameters, i.e., when \t is minimum, for each type of spatial section and matter 
content of the generalized Chaplygin gas model. Ho is given in km/Mpc.s, A in units of c, to in Gy and <Zj 
in units of ao- 

In tables ^ and [21 the values of the parameters for the minimum xl (x 2 divided by the number of SNe 
la) are given for the GCGM with five free parameters (a, Hq, fl m o, ^cOj A) and for other cases where the 
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Table 2: The best-fitting parameters, i.e., when xt is minimum, for each type of spatial section and matter 
content of the generalized Chaplygin gas model using the HST prior. Hq is given in km/Mpc.s, A in units 
of c, to in Gy and a,i in units of ao. 
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Table 3: The best-fitting parameters, i.e., when xt is minimum, for each type of spatial section and matter 
content of the traditional Chaplygin gas model. Ho is given in km/Mpc.s, A in units of c, to in Gy and dj 
in units of clq. 

pressureless mater, the curvature or both are fixed, respectively using the flat prior for H and the HST 
(Hubble Space Telescope) prior for Hq. Analogously, the same estimations are presented in tables and 01 
for the CGM, for up to four free parameters (Ho, Q m o, ^co, A), and in table [5] for the ACDM, for up to three 
free parameters (Ho, fi m o, ^a)- 

It is important to emphasize that, for each case, all free independent parameters are considered simul- 



4 



CGM 
with 




k = 


to — 


n m0 = 0.04 


k = 0, 

^mO = 


k = 0, 
fi m0 = 0.04 


2 


1.1167 


1.1193 


1.1167 


1.1170 


1.1193 


1.1200 


Ho 


65.04 


64.80 


65.04 


65.02 


64.80 


64.76 


^kO 


-0.158 





-0.158 


-0.170 










0.000 


0.000 





0.04 





0.04 


o 


1 1 i^Q 


i nnn 


1 ICQ 
i.lOO 




i 
i 


u.yo 


A 


0.808 


0.813 


0.808 


0.827 


0.813 


0.836 


to 


13.91 


14.02 


13.91 


13.93 


14.02 


14.04 




-0.824 


-0.720 


-0.824 


-0.816 


-0.720 


-0.703 




0.701 


0.697 


0.701 


0.698 


0.697 


0.693 



Table 4: The best- fitting parameters, i.e., when \t is minimum, for each type of spatial section and matter 
content of the traditional Chaplygin gas model using the HST prior. Ho is given in km/Mpc.s, A in units of 
c, to in Gy and a% in units of oq. 

taneously to obtain the minimum of jq. So, for example assuming the GCGM, if we ask for the best 
simultaneous values of (a, Ho, £l m o, ^co, A) then the answer is given by the first column of tabled However, 
in this example, asking for the best value of a by weighing (marginalizing or integrating) all possible values 
of (Ho, fl m o, ^cOj A) yields the estimation in the first column of tabic El whose peak of 0.59 for a is totally 
different from 7.70 as best-fitting parameter ! The parameter estimation issue is addressed by the Bayesian 
statistics of the next section, not by best-fitting in n-dimcnsional parameter space. 

Each column of these best-fitting tables was calculated using a corresponding BETOCS notebook. 
The minimization of xt is obtained in the following way : the initial global minimum taken from the n- 
dimensional discrete parameter space (see next section) is used as initial value to search for the local minimum 
of xt by using the function FindMinimum of the software Mathematica |24j . 

Note that the minimum values for \ 2 using the "gold sample" are worse (i.e., higher, from 1.11 to 1.20) 
than the corresponding ones (between 0.74 and 0.77) using the restricted sample of 26 supernovae |211 122| 
(which have excellent quality). While the minima for \ 2 are always slightly higher when the HST prior is 
used (because the HST prior peak is far from the best-fitting Ho), yielding best-fitting parameter values with 
minor changes. Other important results : the parameter a is usually much bigger than 1, A is often near 
unity, fifco is not far from suggesting a flat spatial section (except for ACDM), f2 m o being null recovers the 
quartessence [27] scenario (except for ACDM) , go and a% point to an accelerating Universe today with 
an age to of approximately 14 Gy. But these results must be compared with a more complete statistical 
analysis to be presented below. 

3 Parameter estimations using BETOCS 

Following the previous works [2011211122] . the Bayesian statistical analysis is employed here instead of the more 
usual frequentist (or standard or traditional) statistics. The Bayesian statistics emphasizes considering only 
the (observational) data you have, rather than simulating an infinite space of data, which is an advantage. 
On the other hand, the Bayesian marginalization process is computationally time-consuming if the number 
of parameters of the theoretical model is large. For the case here, with a maximum of five free parameters 
and low number of data points (157 SNe la), the Bayesian approach is better suited than the frequentist 
statistics. See Refs. (23 12U E01 IS] f° r discussions about the frequentist versus Bayesian statistics and some 
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Table 5: The best-fitting parameters, i.e., when \t is minimum, for each type of spatial section and matter 
content of the ACDM model using flat and HST priors. Ho is given in km/Mpc.s, to in Gy and a* in units 
of do- 

applications in physics. 

The probability of the set of distance moduli po conditional on the values of a set of parameters {pi} is 
given by the Gaussian : 

Pim\{Pi}) « exp^ - ■ ( 4 ) 

This probability distribution must be normalized. Evidently, when, for a set of values of the parameters, the 
X 2 is minimum the probability is maximum. This is a valuable information but is not enough to constraint 
the parameters. 

From the probability distribution a joint probability distribution for any subset of parameters can 
be obtained by integrating (marginalizing) on the remaining parameters, see refs. |21l I22|. So, in order 
to properly estimate a single parameter, the probability distribution must be marginalized on all other 
parameters, usually yielding a quite different result if we try to estimate the parameter in a two or three- 
dimensional parameter space. The reason is that, in such multidimensional parameter space, if a parameter 
has a large probability density but in a narrow region, the total contribution of this region may be quite 
small compared to other large regions which have small probability: in the marginalization process, this kind 
of high PDF region contributes little to the estimation of a given parameter. 

Hence the estimation of a given parameter will be made by marginalizing on all other ones. A detailed 
Bayesian analysis of the independent and dependent parameters is shown in tables Ifil 151 and ITfJI for GCGM, 
CGM and ACDM with flat prior for H , and in tables[7||iand[TT]fbr GCGM, CGM and ACDM with the HST 
(Hubble Space Telescope) prior for Ho- Each column of these tables was calculated using a corresponding 
BETOCS JSj notebook (including the best-fitting calculations for the specific model). 

The following estimation analyses will focus on the tables IBTl 1 II and the accompanying figures. 

3.1 Estimation of H : Hubble is not humble 

The predicted value of the Hubble constant today Ho is the most robust one, with minor changes for the 
different models (GCGM, CGM and ADCM) and cases of fixed parameters. When comparing with ref. |20| . 
Ho is now slightly smaller. 
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Fi gure 1: The one-dimensional PDF for Hq. The thin line shows the HST (Hubble Space Telescope) prior, with the lcr 
(68.27%), 2a (95.45%) and 3<r (99.73%) regions delimited by red, blue and green lines, respectively. The thick line shows a 
typical Ho estimation from SNe la analysis in this article, clearly showing that the HST prior has much larger dispersion and 
its 2(7 region includes the Hq estimation from SNe la. 



Figure ^ shows a typical H estimation from SNe la analysis in this article, and the HST (Hubble Space 
Telescope) prior [25] for Hq, which has a much larger dispersion. The effect of the HST prior on the Ho 
estimation is small : it slightly changes the shape PDF for Ho, the PDF peak moves increases (moves to the 
right) and there are some very small changes on the left and right dispersions. 

The minor effect of the HST prior can also be verified by comparing tables Rjl 151 and ITU1 for GCGM, CGM 
and ACDM with flat prior for H with tables and [fll for GCGM, CGM and ACDM with the HST 
(Hubble Space Telescope) prior for H . 

It is important to emphasize that fixing Ho is not acceptable as the n-dimensional PDF quite depends 
on the Hq parameter. For ACDM as an example : xl lS verv high (1.563), flko = — 0.890l ;24i aim many 
other totally different parameter estimations. 



3.2 Estimation of a 

With respect to ref. |2Jj], the peak values of a are slightly increased and the dispersion is also a little larger, 
for example, the GCGM now gives a = — 0.59io~4i- Imposing that the space is flat or fixing the pressureless 
matter lead to positive best values for a. For example, the quartessence ^3G3 scenario (f2 m = 0) predicts 
a = 0.90j£j§. 

Note that the marginalized estimations differs substantially from those extracted from the minimization 
of x 2 , which gives a large positive best value a, but the dispersion is quite high, so even large positive values 
are not excluded, at least at 2<r level. And figure for the joint probabilities for a and A and figures 01 
and ^ for a and A, respectively, clearly show that the marginalization process changes the peak values and 
credible regions depending on the number of dimensions. 

Of course, the CGM is obtained when a is fixed to unity. From the analysis of the GCGM it can be 
inferred that p(a = 1) = 40.81%, i.e., the CGM is favoured with a probability of 40.81%. Restricting to null 
curvature or fixing the pressureless matter density increases considerably this value, from 68.22% to 91.09%. 
Analogously, the probability to have a > (with more physical meaning) is 66.74% and this value is quite 
increased when one or two parameters are fixed. Both p(a = 1) and p(a > 0) are increased with respect to 
ref. 
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Table 6: The estimated parameters for the generalized Chaplygin gas model (GCGM) and some specific 
cases of spatial section and matter content. We use the Bayesian analysis to obtain the peak of the one- 
dimensional marginal probability and the 2 a credible region for each parameter. Ho is given in km/Mpc.s, 
A in units of c, to in Gy and in units of ao . 

3.3 Estimation of A 

Like ref. |2JJ], the results indicate that the value of A is close to unity, but now the dispersion is slightly 
smaller. In the case of GCGM with no fixed parameters, the marginalization of the remaining four other 
parameters leads to A — l.OOOig'^g. Again, this could suggest the conclusion that ACDM (A = 1) model 
is favoured. However, the accuracy of the computation, due to the step (between 0.01 and 0.02) used in 
the evaluation of the parameter, does not allow this conclusion. Instead, it means the peak happens for 
0.98 < A ^ 1. In fact, fixing the curvature or the pressureless matter, the preferred value differs slightly 
from unity, for example the quartessence scenario, fi m = 0, yields A = 0.987lo'3888- But, differently from 
ref. (20], the CGM now predicts a best value for A smaller than unity, A = 0.860i"o'ogg, and the best value 
becomes smaller when one or two parameters are fixed. 

In figure 0] the PDF for A is displayed, both for the GCGM and the CGM, where the marginalization 
is made in all other parameters. Note that the probability to have I ^ 1 (meaning how much the ACDM 
is rule out) is zero only for the GCGM (due to the step used in the evaluation of A), but this probability 
varies from about 60% to 100% for other cases. 

In figure|21the joint probabilities for a and A are displayed, with a non-Gaussian shape. Comparing with 
ref. | 20| . the peak values now happen for large values of a and A. This figure, compared to figures |3 and 
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Table 7: The estimated parameters using the HST prior for the generalized Chaplygin gas model (GCGM) 
and some specific cases of spatial section and matter content. We use the Bayesian analysis to obtain the 
peak of the one-dimensional marginal probability and the 2 a credible region for each parameter. Ho is given 
in km/Mpc.s, A in units of c, to in Gy and in units of do- 

^ is an illustration of the importance of the marginalization process because it changes the peak values and 
credible regions depending on whether two or one-dimensional parameter space is used. 

3.4 Estimation of f2 m0 and f2 c0 

The unified scenario of quartessence, with no pressureless matter, is again favoured as, for example, the 
GCGM and CGM (without fixed parameters) predict n mQ = 0.000i£ooo and n mQ = O.OOOi^^, respectively. 
The same estimations from ref. [55], Q m o = O-OOIq'qo, and ref. fi m o = O.OOIq'qo, show that the increased 
number of SNe la has substantially decreased the estimated error. Compared to ref. j2D], the dispersion is 
also highly decreased, which once more favours the quartessence scenario. The case of flat Universe has an 
even smaller dispersion for the quartessence model. See figures and |H1 

Like ref. [201; the GCGM, the CGM and the ACDM remarkably lead to almost the same predictions 
concerning the dark energy component, O c o, when all parameters are free: 1.012lQ 48 g, 1.030±g ; 4 g 7 and 
0.976iQ'4 2 g, respectively. By comparing with ref. [20], now O c o (which behaves as fl\ for ACDM) has lower 
values and narrower dispersions, like fl m o, see figures [7HH1 

The joint probability for f2 m o and Sl c o is now smoother for the GCGM and the CGM cases (figure 0, 
with the lcr, 2<7 and 3cr contours of figures and [B] showing significantly smaller regions. So the analysis for 
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Table 8: The estimated parameters for the traditional Chaplygin gas model (CGM) and some specific cases 
of spatial section and matter content. We use the Bayesian analysis to obtain the peak of the one-dimensional 
marginal probability and the 2 a credible region for each parameter. Hq is given in km/Mpc.s, A in units of 
c, to in Gy and ai in units of ao ■ 

ACDM is now quite in agreement with the results of ref. |26|. 

3.5 Estimation of tt k0 

In comparison with ref. |20| , a closed Universe is still clearly favoured, but with slightly smaller probability, 
i.e., p(Qko < 0). But more important, the dispersion for O, k o is now substantially narrowed, as shown by 
figuresEHill The probability to have a flat Universe, p(fl k0 = 0), is now greater (45.70 %, 39.84 % and 15.44 % 
for GCGM, CGM and ACDM), and after setting the pressureless matter it increases, exception being the 
ACDM case. 

With respect to refs. |21l I22j (using the selected 26 SNe la data), the dispersion has also significantly 
decreased, for example : fl k0 = -0.251^;^| and tt k0 = -0.228±g;|j^ versus {l k0 = — 0.74±^ J| [221 and 
n k0 = -0.841^3 [ZD, for respectively GCGM and CGM. 

3.6 Estimation of the age of the Universe, t; 

Due to the larger number of SNe la used here with respect to refs. [21 123- the dispersions have been quite 
decreased. For example, t = 14.42+^ 77 Gy and t = 13.931°;^ Gy estimated here for the GCGM and CGM 
versus t = 15.3+ 4 jG?/ [221 and t = U.2tl 8 5 Gy 22, respectively. 

The predicted age of the Universe when no parameter is fixed has increased with respect to ref. |2U|. 
fortunately not anymore dangerously near the recent estimations age of the globular clusters |32j . to = 
1 12.fi - j 1 
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Table 9: The estimated parameters using the HST prior for the traditional Chaplygin gas model (CGM) and 
some specific cases of spatial section and matter content. We use the Bayesian analysis to obtain the peak 
of the one-dimensional marginal probability and the 2 a credible region for each parameter. Hq is given in 
km/Mpc.s, A in units of c, to in Gy and a% in units of ao- 

3.7 Estimation of the deceleration parameter q 

The values for the deceleration parameter qo are increased (less negative) with respect to ref. [3D]. The 
estimated errors are significantly smaller than the ones of refs. for example : qo = — 0-730^0 328 ano - 

qo = — 0.769±g:|gl versus q = 0.80±g;|| gSj and q a = 0.98±£;g| EH, for respectively GCGM and CGM. 

In all cases, p(qo < 0), the probability to have an accelerating Universe today, is equal to or very near 
100%. 



3.8 Estimation of the scale factor a { the Universe begins to accelerate from 

Another useful quantity is the scale factor at the moment the Universe begins to accelerate, <Zj, keeping in 
mind that the scale factor is normalized with its present value a equal to unity. With respect to ref. |2"U], 
a, decreases when no parameter is fixed. As expected, the larger number of supernovae in comparison with 
ref. [22] yields smaller credible intervals for o», for example a; = 0.626^o' 12 3 a o estimated here for GCGM 
versus a, = 0.67t°;^a of ref. (23- 

The probability the Universe begins to accelerate before today, p{ai < 1) in tables 0-^1 is essentially 
100%, being approximately the same value of the probability to have an accelerating Universe today, i.e., 
p(qo < 0). Theoretically they should be the same, so the fact that these independent probability calculations 
agree almost exactly shows the accuracy and reliability of the Bayesian probability analyses of this work. 
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Table 10: The estimated parameters for the ACDM model and some specific cases of spatial section and 
matter content. We use the Bayesian analysis to obtain the peak of the one-dimensional marginal probability 
and the 2 a credible region for each parameter. H is given in km/Mpc.s, t in Gy and a, in units of oq. 



PDF for GCGM PDF for GCGM with Q m0 = 0.04 




Figure 2: The graphics of the joint PDF as function of (a, A) for the generalized Chaplygin gas model. The joint PDF peak 
is shown by the large dot, the credible regions of 1 cr (68, 27%) by the red dotted line, the 2 cr (95, 45%) in blue dashed line and 
the 3 cr (99, 73%) in green dashed-dotted line. The cases for Q m o = are not shown here because they are similar to the ones 
with Sl m o = 0.04. The cases for fc = are shown in figure 2 of ref. 1201 . 
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Table 11: The estimated parameters using the HST prior for the ACDM model and some specific cases of 
spatial section and matter content. We use the Bayesian analysis to obtain the peak of the one-dimensional 
marginal probability and the 2 a credible region for each parameter. Hq is given in km/Mpc.s, to in Gy and 
ai in units of ao- 

PDF for GCGM PDF for GCGM with f2 m0 = 0.04 




Figure 3: The PDF of a for the generalized Chaplygin gas model. The solid lines are the PDF, the la (68.27%) regions are 
delimited by red dotted lines and the 2a (95.45%) credible regions are given by blue dashed lines. The cases for Q m o = are 
not shown here because they are similar to the ones with Q m o = 0.04. The cases for k = are shown in figure 3 of ref. 1201 . 



4 Conclusions 

The present work has performed the most extensive analysis of the GCGM and CGM in what concerns 
the comparison of theoretical predictions with the type la supernovae data, using the 157 data of the "gold 
sample". By using the high productivity of BETOCS (BayEsian Tools for Observational Cosmology using 
SNe la it was possible to make all the best-fittings, parameter estimations and figures shown here. 

All positive features of GCGM and CGM were enhanced with respect to ref. [201 : 

• the probability to have a > (with more physical meaning) is increased, as well as both p(a = 1) and 
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Figure 4: The one-dimensional PDF of A for the generalized and traditional Chaplygin gas model. The solid lines are the 
PDF, the Ict (68.27%) regions are delimited by red dotted lines, the 2cr (95.45%) credible regions are given by blue dashed lines 
and the 3<r (99.73%) regions are delimited by green dashed-dotted lines. The cases for f2 m o = are not shown here because 
they are similar to the ones with f2 m o = 0.04. The cases for k = are shown in figures 4 and 5 of ref. 1201 . 

p(a > 0); 

• value of A now has slightly smaller dispersion; 

• the dispersion of £l m o and f2 c o are also highly decreased, favouring even more the quartessence |13II27) 
scenario; 

• the dispersion for flko is now substantially narrowed, and the probability to have a flat Universe is now 
greater; 

• predicted age of the Universe (when no parameter is fixed) has increased, fortunately not anymore 
dangerously near the recent estimations age of the globular clusters. 

One important result concerns the Hubble parameter, Hq. It is not acceptable to fix its value because 
arbitrary parameter estimations and usually bad best-fittings are obtained. As the HST (Hubble Space 
Telescope) prior for Hq implies minor effects on all results, we can choose a flat or HST prior. 

The CGM (traditional Chaplygin gas model), where a = 1, remains competitive and preferred in many 
cases : when the five parameters are considered, the probability is 40.81%, but increases as much as to 
95.94% for the quartessence scenario. 

For the parameter A, both GCGM and CGM cases of fixed curvature and matter densities shows a value 
near but small than 1 as the best value of A, such that the ACDM case (A — 1) is almost ruled out. 



14 



PDF for GCGM 



PDF for CGM 




0.1 0.2 0.3 0.4 0.5 0.6 0.1 0.2 0.3 0.4 0.5 0.6 0.7 

!5 m „ 



Figure 5: The graphics of the joint PDF as function of (r2 m o, S7co) for the generalized Chaplygin gas model, GCGM (traditional 
Chaplygin gas model, CGM), where p(f2 m o, £l c o) is a integral of p(ct, Ho, O m o, f2 c o, A) (p(Ho, f2 m o, Q c q,A)) over the (a, Ho, A) 
((Ho, A)) parameter space. The joint normalized PDF peak has the value 4.624 (5.158) for (Q m o,Q c o) = (0.000,0.943) 
((fimOj^eo) = (0.000,1.099)) shown by the large dot, the credible regions of 1 a (68,27%, shown in red dotted line), 2<r 
(95,45%, in blue dashed line) and 3cr (99, 73%, in green dashed-dotted line) have PDF levels of 1.577, 0.345 and 0.038 (2.191, 
0.439 and 0.040), respectively. As f^o + f2 m o + f2 c o = 1, the probability for a spatially flat Universe is on the line Q m o + £l c o = 1; 
above it we have the region for a closed Universe (k > 0, QkO < 0)> an d below, the region for an open Universe (k < 0, fl^o > 0). 



The results indicate that, for the GCGM, CGM and ACDM, a closed Universe is favoured. The GCGM 
and CGM favour the unified scenario (quartessence) , where the pressureless matter density is essentially 
zero. 

There are many current and future applications |201 EH H21 HUD an d developments of BETOCS : using 
different SNe la data sets (SNLS etc), different priors, different cosmological models |55). etc. We 
plan to release variant versions for other observational cosmological data |36l 1371 1551 139| . like BETOCX 
(BayEsian Tools for Observational Cosmology using X-ray gas mass fraction of galaxy clusters), so it will 
be possible to cross the estimations from different observational data. 
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Figure 9: The one-dimensional PDF of Q k0 , Q m o and Q c0 for the ACDM model. The solid lines are the PDF, the la (68.27%) 
regions are delimited by red dotted lines, the 2a (95.45%) credible regions are given by blue dashed lines and the 3a (99.73%) 
regions are delimited by green dashed-dotted lines. The case for f2 m o = is not shown here because it is similar to the one with 
QmO = 0.04. The case for k = is shown in figure If of ref. 1201 . As fi c o = f — £lkO — ^fcO: f° r ^mO = we have f2 c o = f — f2fe0i 
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